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Abstract 



We calculated the differential cross sections of the processes in which one of the pair 
created r particles at an e + e~ collider decays into lepton flavor violating final states e.g. 
r — > fij, t — > 3fi, t — > fiee. Using the correlations between angular distributions of 
both sides of r decays, we can obtain information on parity and CP violations of lepton 
flavor non-conserving interactions. The formulae derived here are useful in distinguishing 
different models, since each model of physics beyond the standard model predicts different 
angular correlations. We also calculate angular distributions of the major background 
process to r — > Z7 search, namely r — ► IvPj, and discuss usefulness of the angular 
correlation for background suppression. 
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1 Introduction 



Recent results form neutrino experiments such as the Super-Kamiokande experiment 
strongly suggest neutrino oscillation so that there are flavor mixings in the lepton sec- 
tor U]. This implies that the charged lepton flavor violating processes such as fi — > ej, 
r — > fi'-f, r — > e7, etc. also occur at some level. It is, therefore, important to search for 
lepton flavor violation (LFV) in rare decay processes of muon and r. 

Prediction of branching fraction of LFV processes depends on models of physics be- 
yond the standard model. In the minimal extension of the standard model which takes 
into account neutrino oscillations by seesaw mechanism of neutrino mass generation 0, 
the expected branching fraction is too small to be observable in near future §. On the 
other hand, in supersymmetric (SUSY) models, the prediction can be close to the current 
experimental upper bound. In this case, the flavor mixing in the slepton mass matrix 
becomes a new source of LFV. Even in the minimal supergravity scenario [|j], in which 
the slepton mass matrix is proportional to a unit matrix at the Planck scale, the renor- 
malization effects due to LFV interactions can induce sizable slepton mixings [||]. For 
example, such LFV Yukawa interactions exist in SUSY Grand Unified Theory (GUT) Q, 
SUSY model with right-handed neutrinos 0, and SUSY models with exotic vectorlike 
leptons |J. Another interesting possibility is models with extra dimensions, where the 
neutrino masses and mixings are obtained from the Yukawa interaction between the ordi- 
nary left-handed leptons and the gauge-singlet neutrinos which propagate in the bulk of 
extra dimensions 0. This Yukawa interaction breaks the lepton flavor conservation and 
the Kaluza-Klein modes of the bulk neutrinos can enhance [i — > decay, r — >• /ry decay, 
etc. through the loop diagrams [|l(J . 

In the muon decay, the polarized muon experiments provide useful information on 
the nature of LFV interactions [pT| , We can define a parity (P) odd asymmetry 

for fi — > ej process and P and time reversal (T) odd asymmetries for fi — > 3e process. 
These asymmetries are useful to distinguish different models. For example, in the SU(5) 
SUSY GUT model with small and intermediate value of tan j3 (a ratio of two vacuum 
expectation values of Higgs fields), only ji + — > e^7 (or jj," — > e^ff) occurs because LFV is 
induced through the right-handed slepton sector. On the other hand, the SUSY models 
with right-handed neutrinos predict /i + — > e^7 (or ^ — > e^7), and the SUSY models 
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with vectorlike leptons can induce both /x + — > e\^ and fi + — > e^7 depending on how 
interaction breaks the lepton flavor conservation. In the models with extra dimensions, 
only fi + — > ejj7 can occur. As for the T odd asymmetry in the fi — > 3e process, it was 



shown that asymmetry could be sizable in the SU(5) SUSY GUT [11], [12 1. 

In this paper, we discuss the LFV processes of r decays such as r — ► /ry, r — ► 3//, 
r — > fiee, taking into account P and T odd asymmetries. In the r + r~ pair production at 
e + e~ collision, we can extract information on the spin of the decaying r particle from the 
angular distribution of the r decay products in the opposite side. Using this technique, we 
can obtain the P and T odd asymmetry defined in the rest frame of r. The method of the 
spin correlation has been developed since the days before the discovery of r particle |TB| . 
There have been many works on spin correlation method in search for anomalous coupling 



involving r [14]]. We have applied the formalism in order to obtain the information on 
LFV interactions under P and T symmetries. We also calculate angular correlation of the 
process where one of the r's decays through r — > Zz/z/7 mode. This mode is a background 
process to the r — > £7 search if the neutrinos carry out little energy. As in the muon 
case |15| ) the angular correlation is useful to identify the background process and the 



background suppression is effective for r~ — > (t + — > fi~Rl) search. 

This paper is organized as follows. In section 2, we introduce a formalism to calculate 
the spin correlation. In section 3, we present a differential cross section of the production 
and decays of r + r~ at e + e~ colliders where one of r's decays in r — > /ry or r — >• ey 
modes, and show how to extract P-odd asymmetry of r decay. In section 4, P and T 
odd asymmetries in three body LFV decays (r — > 3/i, r — > nee, etc. ) are considered. 
In section 5, we consider r — ► Iv&j mode and show that the analysis of the angular 
distribution is useful for background suppression of the r — > fij and r — >• ey searches. A 
summary is given in section 6. Appendices contain the derivation of basic formulae and 
a list of kinematical functions. 



2 General formula for spin correlation 

In this section, we present general formulae used in the calculation of differential cross 
sections and spin correlations. 

We calculate differential cross sections of e + e~ — > t + t~ — > /b/a, where fs (Ja) repre- 
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sents the decay products of r + (t~). If the intermediate states were spinless particles, the 
cross section is simply a product of a production cross section and decay branching ratios. 
However, in the case of spin 1/2 particles, we have to take into account spin correlation 
between two intermediate particles. If we take r + — > f B to be a LFV decay mode, we 
can measure P and T violation of LFV interactions by using the angular correlations of 
decay products of r + and r~ . 

The differential cross section of e + e~ — > t + t~ — > JbJa is given by 

3 

da = da p dB T ~^ fA dB T+ ^ fB + dI] ab dRl~ tfA dR T b + ^ fB , (1) 

a, 6=1 

and 

p d 3 pA d 3 pB 1 /„ u j4/ x p / n x 

da = j^m j^m 2s (2?r) 6 <pa+p*-p*-p<-) a > ( 2 ) 



iTDT-^fA 1 < ^ 3 ^ 1 ^ 3 ^ n * Co x4 ( \ D- 



T (2vr)32g0 +1 (27r)32g° +m 2m 



(27r) 3 2p0 (2vr) 3 2^ 2s 



(3) 



1 / ra+m \ 

-(2vr) 4 5 4 ^ a D +, (4) 

™ T \i=n+l / 



^ = To^T To^flT o 1 ( 2 ^) 4 *V +PB- Pe+ - p.-) P P , (5) 



(6) 



dR « =f wm'''wm^ { ] 5 \u qt ' PA ) pa ' 

where we assume that is a n body system and f B is a m body system. p e + (p e -) is e+ 
(e~) four momentum, p B (j>a) is r + (r~) four momentum, and ^'s are momenta of final 
state particles, s is determined as s = (p e + +p e -) 2 - T and m T are the width and the mass 
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of the r, respectively. In order to define a p , a D ~ , a D+ , p p b , p®~, and p b + , we first write 
down the invariant amplitude of e + e~ — > t + t~ — > /b/a as follows: 

e 2 - 

Af = — A(p A + m T )^(p B ~ m T )B v e +7 M u e - 
11 

x ■ 



Pa~ ( m r~ f) 2 Pi - K- f) 2 ' 
where u e + (it e - ) is the wave function of positron (electron) and A and B are spinors which 
include wave functions of final states and interaction vertices. By using the Bouchiat- 



Michel formulae [16] and the narrow width approximation, a , a , a + , p p b , p^~ , and 
are given by 

1 e 4 

a p = - - Tr + m T )^{p B - m T )Y] Tr &+7/<A»-7*] , (9) 

cP- = l - {A{p A + m T )A} , (10) 

a D+ = \ {B(p B -m T )B} , (11) 

p P ab = \j 2 Tr [75 & (Pa + m T )7^75 A (1>b - m T )Y] Tr [& +7 ^ e -7,] , (12) 
P D a - = \ {A l5 &(p A + m T )A} , (13) 

P° + = \ {B {j>B - mr)B) , (14) 

where the spins of the final state fermions are summed over, and four vectors {s a A )^ and 
{s B ) u (a,b = 1,2,3) are a set of vectors which satisfy following equations. 

PA-s A =p B -s b B = 0, (15) 

s A -s\ = s B - 4 = -S ab , (16) 

5>aW4), = + , £(4U4), = -g„ + VJ ^ • (17) 

a=l r 6=1 r 

The derivation of the above result is shown in Appendix [X]. Notice that da p , <i-B T 
and dB T+ ^f B in eq.(|I]) are the t + t~ production cross section and r decay branching 



ratios, in which spins of r's are averaged, and dH^ b , dR T a ^^ A , and dR T b s represent the 
spin correlation effects of this process. 

In above formulae, it is assumed that the r pair production occurs through the photon 
exchange. It is straightforward to include the contribution from the Z boson exchange 
and the 7 — Z interference. If we consider the e + e~ center of mass energy to be in 
the range of the t + t~ threshold energy considered in the r-charm factory or the T(4S) 
resonance energy where e + e~ B-factories are operated, these effects only contribute to the 
production cross section at the level of O(10~ 4 ) of the photon exchanging diagram. 

3 Parity asymmetry in r — > fij decay 

Let us calculate the cross section of e + e~ — > t + t~ — > /i + 7 + f& processes. For /a, we 
consider hadronic and leptonic modes such as (tcu, pv, a\v, and Ivv). Below we neglect 
the muon mass compared to r mass, and therefore all formulae can be applied also to the 
t — ► ey process. The effective Lagrangian for r + — ► p +/ y decay is given by 

C= /| {m T A R fa^P L pF^ + m T A L fa^P R pF^ + h.c.} , (18) 

where Gf is the Fermi coupling constant, Pl — (1 — 7s)/2, and Pr — (1 + 7s)/2. The 
operator with the coupling constant Ar (Al) induces the r + — > ji]ff (r + — > 7x^7) decay. 
As mentioned in Introduction, each model of the physics beyond the standard model pre- 
dicts a different ratio of Al and Ar. For example, the SU(5) SUSY GUT in the minimal 
supergravity scenario predicts that only Al has a non- vanishing value for small and inter- 
mediate values of tan/3. Therefore the separate determination of Al and Ar provides us 
important information on the origin of LFV. For this purpose, we need information about 
the t polarization. This can be done by observing angular distributions of final state of 
r decay in the opposite side in the modes of r — > m/, t — > pi/, t — > a\v, and r — > Ivv, 
because these processes proceed due to the V — A interaction and therefore have a specific 
angular distribution with respect to polarization of r. Using r + — t~ spin correlation, we 
can determine \Al\ 2 and |Ar| 2 , separately. 

We first define three coordinate systems (Fig.|I|). The first coordinate system (Frame 
1) is the center of mass frame of the e + e~ collision in which the z axis is taken to be the 
e + momentum direction. The second one (Frame 2) is the rest frame of the r + , and the 
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Frame 2 




Ui Frame 1 



Figure 1: The coordinate systems. The plane determined by e + e and r + r momentum 
vectors corresponds to the xz-planes in each of three coordinate systems. 

third one (Frame 3) is the rest frame of the r~ . More explicitly, the relation of a four 
vector in the three systems are given as follows: 
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& ,(19) 



where 7 = v / s/(2m r ) and (3 T = ^Jl — 4m^/s, and the four vectors £i_3 are defined in 
Frame 1-3, respectively. We calculate the production process in Frame 1, and the r + (t~) 
decay in Frame 2 (Frame 3). In the calculations, we choose the spin vectors (s^) M , {s^)^ 
as follows: 



[s a A r 







5, 



a+1 



[in Frame 3) 



(20) 



:4r 







(in Frame 2) . (21) 
The production cross section and spin dependence term are obtained from eq.@ and 
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eg. fll2|) as follows: 



do 1 



dQr na 2 



dZ 



ab 



An s 
dQ T na 2 

An S 



Ami 



1 + 



Ami 



Ami i 

T - cos 2 6 T 



(22) 



Ami 



( (l + ^)sin 2 # T 



x 



^sin20 T 





s 







1 _l 4m| 



sin 2 6> 



1 _ 



cos 2 9 T 



(23) 



where # T is the angle between the e + and r + directions in the Frame 1, and dQ T is a solid 
angle element of r + , <if2 r = dcos6 T d(j) r . 

For decay processes, we take r + — > /i + 7 for the r + side and hadronic (r~ — > n~v, 
t~ — > p~v, and r~ — > aiu) and leptonic (r~ — > l~vv) decays for the r~ side. dB T+ ^^ +1 
and dR[ 7 (see eq.(P) can be calculated from eq. flTTD and eq. flU]) in which the spinor 



B is given by 



B = -= G F m r0 r^(g 7 ) M (A«P L + A L P R )elv{q^ 



(24) 



where is the polarization vector of the photon and (g 7 )^ is the momentum of the photon 
and v (g M ) is the wave function of the muon. These quantities are given as follows: 
dQ„ 1 2 



dB 1 



^ T 7 



An T Ti 



Glml (\A L \ 2 +\A R \ 2 ) 



(25) 



" ' " Glml {\A L \ 2 - \A R \ 2 ) 



An T n 



sin 6^ cos M 
sin ^ sin M 
cos#„ 



(26) 



where 



A" VAt, 



are angles in the polar coordinate for a unit vector of the muon momentum 



direction in Frame 2. The three components in eq.(p6|) corresponds to b = 1, 2, and 3. 
Next, we list dB and dR for the r~ decay in each mode of t~ — > 7r~z/, t~ — > p~i/, 



r — > a x i/, and r — > / w. For r — ► 7r i/ decay, the spinor A in eqs.(|T0|) and (p~3[ ) is 
given by 



A = 2iV ud f 7r G F fa w P L u(q^ 



(27) 
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where f n is the pion decay constant, q n is the momentum of the pion, and u(q u ) is the 
neutrino wave function. Then, dB T v and dR T a v are given by 

dir--*-" = ^ i -L |K,| 2 /. 2 G|m 3 r , (28) 

(sin 6^ cos 0„- \ 
sin ^ sin 0^ , (29) 
cos 9 W J 

where (6 n , <j) n ) are the polar angles of tt~ momentum in Frame 3 and dQ n = dcos8 n d^^ 
(we use a similar notation in the following expressions). Here we neglect the mass of the 
pion. As before three elements in eq. (|29|) corresponds to a = 1,2, and 3. Similar results 
can be obtained for the vector mesons. The spinor A for r — * pu, r — > a\v is given by 

A = -2V ud g v G^ v P L u(q v ) . (30) 

where gy and ey are the decay constant and polarization vector of the corresponding 
vector mesons, respectively. From this expression, we can obtain dB and dR for the 
longitudinally polarized vector mesons e.g. r~ — > p~[V)v and r~ — > a{{L)v as follows: 

dBT^r* = ^lil \V ud \ 2 ( ^) 2 GW T mUl ~ 4) , (31) 



4/T 1 671 \ m V r / m T 

(Sin #y COS 0y \ 
sin#ysin0y , (32) 
cos 9 V J 

where my and (#y, <j>v) are the mass, and polar angles of the corresponding vector meson, 
respectively. For the transversely polarized vector mesons, the spin dependence terms 
have a minus sign contrary to the case of the pion and longitudinally polarized vector 
mesons. 

tr-wr. = *?v i _L „y> (9l\' gWMAi _ , (33) 

47T 1 87T V 77 ^ / m T m T 

sin #y cos 4>y 

dR T ~^ v{T) ~ v = dB T ~^ v{T) ~ u \ -sin^ysin0y ] . (34) 

COS 9y 
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For leptonic decays, after integrating over the phase space of the neutrinos, the branching 
ratio and the spin dependence term are given by 



dRl 



->/ vv 



->/ vv 



dn t , 1 G\m\ o , 
ax — — 



An 

An 



r !92tt 3 



2x 2 (3 - 2x) 



1 r*2„,5 / sin^cos^ 

dx 1 Gp ^ -~ 2 ' 



r i92tt 3 



I 2x 2 (l - 2x) I sin 6i sin . 

COS 



(35) 
(36) 



where we neglect the masses of the leptons (e and fi) and x is the lepton energy normalized 
by the maximum energy m T /2 i.e. x = 2Ei/m T , and (6i,4>i) are the polar angles of the 
lepton in Frame 3. 

Substituting these results into the formula in eq.(0), we obtain the differential cross 
sections of each processes. For example, the differential cross section of e + e~ — ► t + t~ — > 
/i + 7 + n~v process is given by 



da(e + e — > r + r — > |U + 7 + n u) 
<j(e + e~ — > r + r~ x 



B( 



7T 



i/) 5(r+ 



/^ + 7) 



dQ T dVt^ dQ n 



An An An 



x 



1 + 



Amt 



Ami ' i 
1 T - ) cos 2 9 T 



+Ap ( sin 6^ cos 0^ sin 9 n sin 0^ cos 6^ ) 



2m T 



sin 26 T 



x 



%sin20 T 

sin 6* M cos /t 
x | sin 6^ sin </> M 



1 - 



but 



sin # T 







- 1- 



1 + 



cos 2 9. 



cos 



where 



trie e 



r r 



4m?. / 2m 2 
1 1 1 + - 



3s V s \ s 

is the t + t~ production cross section. The branching ratio of r~ — > n~v and r H 
is given by 

1 1 



B{ 



n 



r 8 



n 



\V ud \ 2 flG 2 F ml , 



(38) 
/i+7 

(39) 
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B{r+ -> ^7) = I ~ G 2 F m 5 T (\A L \ 2 + \A R \ 2 ) , (40) 
and the asymmetry parameter Ap is defined as follows: 

\a l \ 2 + \ArP 

We can see that the measurement of angular correlation of the pion and muon momen- 
tum enables us to determine the parameter Ap, so that we can obtain \Al\ 2 and \Ap\ 2 
separately. 

A simpler expressions can be obtained if we integrate over the angle 9 T , (fi T , (fi n , and 
7 in eq.([37|). The differential cross section is given by 

d<j(e + e~ — > t + t~ — > /i + 7 + n~v) 

t + - + -\ui + + \nt - - ,dcos9 fl dcos9 7r 

= a(e^e — > t^t )B(t^ — > ^T^jByr — > n v) 

2 2 

(s — 2m 2 \ 
l- 7T T i ^Ap cos 9, cos 9 W J . (42) 

Notice that angular distribution in the rest frames of r + and r~ can be easily converted 
to the energy distribution in the center of mass frame of the e + e~ collision. We obtain 

dcr(e + e~ — > t + t~ — > jU + 7 + n~v) 

s 

= <j{e + e~ — > t + t~)B{t + — > fi +/ y)B(r~ — > 7T~v) — - dz^dz^ 

s — Ami 



x ( 1 - ~ ZZ^ 7o„9N a p - !)( 2 ^ - !) ) > ( 43 ) 



s(s — 2m 2 / 
(s — 4mJ)(s + 2m 2 ) 

where z^ = E^/E T (z n = E n /E T ), and E^, E n , and E T = \fsj1 are the energies of the 
muon, pion, and r in the center of mass frame, respectively. 

The angular (or energy) distributions in eq.(^2|) (eq.(f43|)) can be understood as follows. 
Because of the helicity conservation of the t + t~ production process, the helicities of r + 
and t~ are correlated, namely t^t^ or t^t^ is produced. This means that two r spins 
are parallel in the limit of \fs 3> m T . In the decay process, the tt~ tends to be emitted to 
the spin direction of r~ for t~ — > tt~i>, because of the V — A interaction. On the other 
hand, for r + — > /i + 7 decay, the muon tends to be emitted to the same direction of the 
r + spin if Ap > 0. Therefore the differential branching ratio is enhanced (suppressed) if 
the sign of cos 9^ cos 9 n is negative (positive). In other words, pion and muon energies in 
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Figure 2: The weight function w(u,v). 

the center of mass frame of the e + e~ collision have a negative correlation if Ap > 0. If 
Ap < 0, we have an opposite correlation. 

We can define an asymmetry A At+7 ' 7r v by the following asymmetric integrations 

f dcos9 u dcosO^ w(cos6L, cos6L) — ^? — «- 

Ajj, + JtW - v J " ^ d cos O^d cos W 

cr(e + e~ — > t + t~)/3(t + — > fi +/ -f)B(r~ — ► 7r~z/) 
AT++ + TV - iV+~ - /V~+ 
~ TV++ + TV— + 7V+- + 7V-+ ' ^ ^ 

where the weight function w(u,v) is defined by 

, . uv . . 

V J = 1 — r , (45) 

and shown in Fig.§. In the second line, 7V ±=t are the event numbers where the first ± 
represent the sign of cos 9^ and the second one is that of cosO^, respectively. A^ 1 ^ v is 
related to the parameter by 

In Fig.[3], y/s dependence of A^" 1 ^ 7 is shown for Ap = —1. We can see that the asym- 
metry is already close to the maximal value at the B-factory energy. 
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Figure 3: The observable asymmetry y 4 /i+7,7r v vs y/s for Ap = — 1. The dashed line 
represents the \fs of r-charm factory and the dotted line represents that of B-factory. 

It is straightforward to extend the above formula to other cases. We only present here 
formulae corresponding to eq.(^) for different decay modes of r~ . 

da(e + e~ — > t + t~ — > p, +r y + V~v) 
= a(e + e- -> t + t-)B(t + -> ^)B(t~ -> V~u) 



d cos 9^ d cos dy 



S — 2?TJ 2 

x ( 1 ± — — — T - Ap cos 9^ cos 9 V ) , (47) 



s + 2w? T 

where + corresponds to the vector mesons with transverse polarization V = p(T), ai(T) 
and — corresponds to those with longitudinal polarization V = p(L), ai(L). For leptonic 
decay, we obtain 

da(e + e~ — > t + t~ — > /i + 7 + Z _ z/z/) 
= a(e + e- -> t-+t-).B(t + -> / u + 7) J B(r" -> r^) cfe 2x 2 

f s — 2m 2 1 

x <^ 3 - 2x m - 2x)A P cos 0„ cos Q x \ . (48) 

[ s + 2m~ J 

The measurement of the polarization of the vector mesons can be done by the analysis of 
the distribution of the two (or three) pions from the p (ai) meson decay [|17j . 
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In the case of r decays into /x 7 and r + decays via V — A interaction, the dR T a ~*$ A 



+ f 

and dPZ acquire extra minus signs. For example, 



dR r--*»-y = __M _ _ G 2 m 5 (^2 _ ^2) [ sin ^ sin ^ | f (49) 



cos^:, 



= -f ^ iKjV^m* | sin^sin^; | , (50) 



sin 6^ cos d>'„ 
i ^ sin 
cos 9' 



where (9' <f>' ) ((9' n , <f>' n )) are the polar angle of the muon (pion) momentum in Frame 3 
(Frame 2). The formula in eq.(^) can be applied to the r~ — > fi~j case by the replacement 
of (9^,9 n ) by (9' ,9!^), and therefore same angular and energy correlation holds as in the 
t + — > /i + 7 case. In a similar way, we can obtain the formulae corresponds to eqs.(f47|) and 
(|8|) for the r~ — > case by replacement of (9y,9[) by (9' v ,9i), where 9' v {9[) is the 
angle between the vector meson (lepton) momentum and r + direction in Frame 2. 

4 P and T asymmetries in LFV three body r decays 

In this section, we consider LFV three body decays i.e. r — > 3/i, r — > 3e, r — > /zee, and 
t — ► e/z/z. Within the approximation that the muon and electron masses are neglected, 
t — > 3/i and r — > 3e (or r — > /zee and r — > e/i/i) give the same formula, so that we only 
consider r — > 3/i and r — ► /zee processes. In these processes, we can define the P odd as 
well as T odd asymmetries of r decays. 

For r + — > /i + /i + /i^ decay, the effective Lagrangian is given by 

£ = --^{m T A R fa^P L fiF^ + m T A L fa^P R fiF^ 

+9i (fP L fJL) {fiPLfi) + 92 (fP R ^) (fiPnfJ.) 
+^ 3 (r7 M P R /i)( / u7 M P /? /i) + g^ffPL^ifirf^PLli) 
+g 5 (fYP Rf i)(fi^P L fi) + g 6 (f^P L ft)(fi lfl P R fi) + h.c.} . (51) 

With this Lagrangian in eq.([5l|), we can calculate the differential branching ratio <iB r+_+3M 
and the spin dependence term dRl + ^ Sfi in eq. ([!]). In order to calculate these quantities 
we first define the Lorentz frame (Frame 4) for the three body decays fL2 |. Frame 4 is 
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Figure 5: The relation between Frame 2 and Frame 4. 
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the rest frame of r + and we take z-direction to be \f momentum direction, and xz-plane 
to be the decay plane. The x-direction is determined so that the a;-component of the 
momentum for the fi + with larger energy is positive. The coordinate system are shown 
in Fig.f|. Any four vector in Frame 4 is related to that in Frame 2 by the Euler rotation 
with three angles (9, 0, ip) as follows (Fig.||): 



( 1 o \ / 1 

cos — sin 

sin0 cos0 

\ 1 J 







\ / 1 



cos 6* sin# 
10 
\0 -sin6 cos 9 J 







cos-0 
sin-0 
\0 







— sin ip 

COS?/' 

1 / 



S2 > 



(52) 



where 



< < vr , < < 2vr , 



< ip < 2tt . 



(53) 



We also define the energy variables X\ = 2E\/m T and x 2 = 2E 2 /m T where E\ (E 2 ) is the 
energy of /i + with a larger (smaller) energy in the rest frame of r + . 

With these angles and energy variables, the branching ratio and spin dependence term 
can be expressed as follows: 



1 m 5 Gr 2 

^ B r = t — F ^ ^ j cog q 

r 2567T 5 



dip X 



(54) 



1 m 5 G 2 

— - dxi dx 2 dcosO dd) dip 



r 256tt 

-Ys e + Z^CgC^ - S^Sip) + W^CgS^ + C,^) 

x J Ys e s$ + Z{-cgc^s^, - s$c$) + W^-cgs^s^ + c^) 
Yc e + Zsgc^ + Wsgs^ 



(55) 



where sg (s$, s^) and eg (c$, c#) represent sin 9 (sin0, sin-0) and cos6> (cos0, cos-0), 
respectively. The functions X, Y, Z, and W are defined as follows: 



X 



gJ2 1-12 

' 1 + -ttt + \93\ 2 + M 2 ) a 1 (x 1 ,x 2 ) + (\g 5 \ 2 + Ige] 2 ) oc 2 {x 1 ,x 2 ) 



16 16 

+ (leA^I 2 + \eA L \ 2 ) a 3 (xi,x 2 ) - Re(eA R gl + eA L g* 3 ) a 4 (xi,x 2 ) 
-Re(eA R g* + eA L g* 5 ) a s (x 1 ,x 2 ) , 



(56) 
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Y 



\9i\ 



92 , | ,2 



16 16 

-Re(eA i? 5g - eA L g* 5 ) a h {x x ,x 2 ) + (\g 5 \ 2 
+ (\eA R \ 2 - \eA L \ 2 ) p 2 (x 1 ,x 2 ) , 



g 4 \ 2 a 1 {x 1 ,x 2 ) + Re(eA R gl - eA L gl) a 4 (x 1 ,x 2 ) 



2 ) fii(x li x 2 ) 



(57) 



z = {\9b? ~ Ifl'el 2 ) Ji{xi,x 2 ) + (\eA R \ 2 - \eA L \ 2 ) 72(2:1, x 2 ) 

-Re (eAj^| - eA^g) 73(2:1, £2) + Re (eA R g* 6 - eA L g\) 74(2:1, x 2 ) , (5£ 



-Im (eA^ + eA L g* 3 ) 73(2:1, ac 2 ) + Im (eAj#g + e^L^) 74(^1^2) 



(59) 



where e (> 0) is the positron charge and functions cti-5, /3i_ 2 , and 7i_ 4 are given in 
Appendix [FJ. Notice that the Y and Z terms represent P odd quantities with respect 
to the r + spin in the rest frame of r + and the W term represents a T odd quantity. 
These are the same as P and T odd terms considered in the differential decay width of 



e T e T e 



3. 



The differential cross section is obtained by substituting this into eq.(|l]). In the case 
that the opposite side r decays into tt~u, we obtain after integrating over (p n , if), 9 T , and 



da(e + e 
= a(e + e" 



r r 



jj, [I jj, + 71 V) 



t + t-)B(t- 



COS6L , 

k v) I V I dx\ dx 2 dcosO 

128tt 4/ 2 



X — I {Y cos 9 + Z sin 9 cos + W sin 6* sin 0} cos 9 m 



(60) 



s + 2m 2 . 

The terms X, K, Z, and W can be extracted by the following (asymmetric) integrations. 

d 5 a 



dcos9 dd> dcos9 m 



dx±dx 2 d cos 9 d<p dcos9 7T 



x < (j[e e 



t + t-)B(t- 



ir v) 



m\G\ 
32tt 3 



/r 



-1 



x 



(61) 



d cos 6* d cos 9 n w (cos 0, cos 9 n 



dx\dx 2 dcos9 dcos9 7T 
wtG 



x iree 



r r 



£G| /r \r J (.s-2/;r; 



32tt 3 



4(s + 2m2) 



(62) 
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dcosOn w(coscf), cos^) 



dx\dx2 d(j) d cos 



^(^..VMr.oKf/r)}" ="i^§ z - (63) 



<i cos 6L wfsin cos 0„ 



dx\d%2 d(f) dcos9 n 



x { ff(e+e - W (^/r) = . (64) 

Notice that the function W represents CP violating LFV interaction. We can see that 
this is induced by the relative phase between the photon-penguin coupling constants (Al 
and A R ) and the four-fermion coupling constants (#3 — g§). 

A similar formula can be obtained for the r + — > fi + e + e~ decay. The effective La- 
grangian for the r + — > fi + e + e~ is given by 

C = ^{m T A R fa^P L iiF^ + m T A L fa^P R ^ u 

+\ 1 (fP L fi)(eP L e) + \ 2 (fP L fi)(eP R e) 

+X 3 (fP Rf i)(eP L e) + X 4 (fP Rf i)(eP R e) 

+\ 5 (fYP Ll i)(e lfl P L e) + \ 6 (f^P Lf i)(e^P R e) 

+X 7 (f^P Rf J,)(e^P L e) + Xs(f^P R fi)(e^P R e) 

+\ 9 (fa^P L [i)(ea flu e) + \ lQ {fa^ P R ^){ea^e) + h.c.} . (65) 

In this calculation, we define Frame 4' which is almost the same as Frame 4 in the r + — > 

[i + fi + fi~ case. The definition is obtained by the replacement of /i~ of r + — > fi + fi + fi~ by e~ 

of r + — > yU + e + e~, the /i + with a larger energy of r + — > fi + fi + fi~ by fi + of r + — > fi + e + e~, 

and n + with a smaller energy of r + — > by e + of r + — > /i+e+e - . If we take the 

definition of (9, <fi, ip) in such a way that the same relation is satisfied as in eq. (|5^) , the 

branching ratio and the spin dependence term are given by 

dBT +->n+e+e- = 1 d d d q d ± d ^ x , 

T 256tt 5 ^ y K ' 

1 772^ (jr^ 

dRl ^ M e e = 1 — ^ cfei dxi dcosO dd dib 

b T 256vr 5 ^ ^ 

-y's C^, + Z'^CgC^ - S^) + W'tcgS^ + C^) 

Y'sg s 4 , + Z' (-cec^s^ - s^) + W'(-c e s^s^ + c^) 
Y'c e + Z'soCf, + W'sgs^ 

(67) 



18 



where the functions X', Y', Z', and W are given by 

X' = (\eA R \ 2 + \eA L \ 2 )A 1 (x 1 ,x 2 )+Re(eA R X; + eA L X* s )A 2 (x 1 ,x 2 ) 

+Re (eA R \; + eA L X* 7 ) A 3 (x 2 ) + (|A 4 | 2 + \X 2 \ 2 + |A 3 | 2 + |A 4 | 2 ) MM 
+ (|A 5 | 2 + |A 8 | 2 ) M(xi,x 2 ) + (|A 6 | 2 + |A 7 | 2 ) A 6 (x 2 ) 

+ (|A 9 | 2 + |Aio| 2 ) A 7 (x 1: x 2 ) + Re (A1A9 + A 4 A* ) MM x 2 ) , (68) 

Y' = -Re(eA R X* 5 -eA L X* 8 )A 2 (x 1 ,x 2 ) +Re(eA R X* 6 - eA L X* 7 ) A 3 (x 1 ,x 2 ) 
- (|A 5 | 2 - |A 8 | 2 ) A h (x u x 2 ) + (\eA R \ 2 - \eA L \ 2 ) B 1 (x 1 ,x 2 ) 
+ (|Ai| 2 + |A 2 | 2 - |A 3 | 2 - |A 4 | 2 ) B 2 (x u x 2 ) + (|A 6 | 2 - |A 7 | 2 ) B 3 (x h x 2 ) 
+ (|A 9 | 2 - |A 10 | 2 ) B 4 (x u x 2 ) + Re (A1A9 - A 4 A* ) B 5 (x u x 2 ) , (69) 

Z' = (|e^| 2 -|eA L | 2 )C 1 (x 1 ,x 2 ) + Re(eA iJ A;-eA L A*)C 2 (x 1 ,x 2 ) 

+Re (eA R X; - eA L X* 7 ) C 3 ( Xl ,x 2 ) + (lA^ 2 + |A 2 | 2 - |A 3 | 2 - |A 4 | 2 ) C 4 (x h x 2 ) 
+ {|A 6 | 2 - |A 7 | 2 + Re (-2AxAS + 2A 4 A* )} C 5 {x lt x 2 ) 

+ (|A 9 | 2 -|A 10 | 2 )C , 6 (x 1 ,x 2 ) , (70) 

W = lm(eA R Xl + eA L X* s )C 2 (x u x 2 ) + lm(eA R X; + eA L X* 7 )C 3 (x u x 2 ) 

+Im (AiA; + X 4 X* 10 ) C 7 (x u x 2 ) . (71) 

The functions A 4 _ 7 , -Bi_ 4 , Ci„ 6 are given in Appendix [B]. The X', Y', Z', and W can be 
extracted in the same way as in eqs.(|6T])-(|6"4"|). 

Next we consider the decay mode of r + — > fi~e + e + . This case is different from above in 
the point that both r — ► e and \i — > e transitions are necessary. The effective Lagrangian 
for this process is given by 
AG 

£ = --^{g' x {fP L e){flP L e)+g' 2 {fP R e){flP R e) 

+g' 3 {f 1 »P R e){iil»P R e) + g' A (fYP L e)(m»P L e) 
+g' 5 (frP R e)(m,P L e) + ^(^P^)^^) + h.c.} . (72) 

If we take a coordinate system similar to Frame 4, in which the larger (smaller) energy 
fi + is replaced by the larger (smaller) energy e + , dB T+ ^^ e+e+ and dPJ b ^ M e e are given 
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by 

dB r+^-e+e+ = * ^ ^ rf cog £ ^ ^ (73) 

T 2567T t ' 

1 TTl^ 

dRl ^ M e e = — ^ | rfxi cfx 2 dcosO dp dip 

(-Y"s e + Z" {coc^ - s^) \ 
F"s fl + Z"(—C0C^ - s^Cip) , (74) 
y% + z% c ^ / 

where functions X", y", and Z" are given by 

X" = (^ + ^ + |^| 2 + l^| 2 )«i(^i^ 2 ) + (l^| 2 + KI 2 )«2(x 1 ,x 2 )^ (75) 

Y" = (^-^ + l^| 2 -|^| 2 )«i(x 1 ,x 2 ) + (|^| 2 -|^| 2 )A(x 1 ,x 2 ), (76) 

Z" = (l^| 2 -K| 2 )7i(x 1; x 2 ), (77) 

where ai_ 2 , /3i, and 71 are the same functions that we defined in r + — > fi + fi + fi~ calcu- 
lation. X", y", and Z" can be extracted by asymmetric integrations as before, but we 
cannot obtain information on CP violation in this case. 

Notice that the above three cases exhaust all possibilities in the three body decay of 
r to e and/or \i as long as we neglect the electron and muon masses compared to the r 
mass. Namely, the formula for other cases can be obtained by appropriate replacements 
of e and/or fi. 

The formulae for LFV decays with r~ can be obtained in a similar substitution as the 
t — ► /i7 case. Using appropriate angles of r~ decay in Frame 3 and r + decay in Frame 2, 
dRb gets an extra minus sign in eqs.fl5"5]), QB7|), and (|7ip. 



5 r — > /xz/z/7 process and background suppressions 

In this section, we consider the background processes for the r — > /ry search, and we show 
that the measurement of angular distributions is useful in identifying the background 
process. In the muon decay, the physical background can be suppressed if we use polarized 
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muons [15|. In the following, we show a similar suppression mechanism holds for r decay 
if we use the spin correlation. 

One of the main background for the r — > /ry search comes from the kinematical 
endpoint region of the r — > fivv'-f process where two neutrinos carry out a little energy 
at the rest frame of r. In the following, we assume that r + decays into y^vv^i and 
t~ decays through one of hadronic and leptonic decay processes. For the r~ decay, the 
differential branching ratio and the spin dependence term are given in eqs.(|28"|)-(|36|). For 
t + — > /i + z/z/7, these quantities are given by 

dB B.G. = 1 G ? m } a dx dy dz dQ. sinz — F , (78) 

T 3 X 2 li 7T 5 ?/ 

j nB.G. 1 Glm^a AAA ✓vo ^ 

ait,, = — - — - ax ay dz aSL smz — 

b T 3 x 2 n 7r 5 y 

(sin 9 M cos M \ 
sin ^ sin M , (79) 
cos M / 

where x and y are the muon and photon energies normalized by m T /2, respectively, and 
{6 hi 0m) * s the P°l ar coordinate of the unit vector of the muon momentum direction, all 



defined in the rest frame of r + (Frame 2). /3 M = \J\ — Ar/x 2 with r = m^/m 2 .. The 
angle z is defined by z = 7r — 0^, where is the angle between the muon and photon 
momentum in the same frame. These quantities can be obtained by a simple replacement 
from the formula of the differential decay width for the radiative muon decay presented 



in Ref . |1| . For completeness, the functions F, G, and H are given in Appendix ||. 

The background comes from the kinematical region near x = 1 + r and y = 1 — r, 
at which the branching fraction vanishes. However, with finite detector resolutions, this 
kinematical region gives physical backgrounds. If we take the signal region as 1 + r — 5x < 
x < 1 + r and 1 — r — 5y < y < 1 — r, the leading terms of the branching ratio and spin 
dependence term expanded in terms of r, Sx, and Sy, after integrating over z, are given 
by 



dB B. G . „ I Glmla ( 4 2 8 3 3 » 

T 3 x 2 n 7r 5 M V 3 1 



dRf- G - * ± j^ML d Q, I -5*% 2 + | 5x 3 5y s ) ( sin £ M sin ^ | . ( 81 ) 



cos Op 
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Then after integrating over M , <f> n , <f> T , and 9 T , the differential cross section for e + e - 
t + t~ — > /i + z/z^7 + 7c~u is given by 

da(e + e" — > t + t~ — > y^vv^ + 7T~v) 

, , _ — \ f Gpin^a \ d cos 9 ^ dcos9 v 

= (Tie e — > T T )BiT — > 7T z/ — r/1 

v J y 1 \3 x 2% 4/ y 2 2 



(82) 



If the photon energy resolution is worse than the muon energy resolution, the term Sx 4 5y 2 
is small compared to (8/3)5x 3 5y 3 . In such a case, the angular distribution is similar to 
the Ar = 0, Al 7^ case of the r — > fij angular distribution. See eqs.(E5|), fl2"6|), and 
(|2|). This feature is useful for the background suppressions for r + — > Hjff search because 
signal and background processes have different angular correlation. For r + — ► ^7 search, 
the signal to background ratio is almost the same even if we take into account angular 
correlation. 

A similar background suppression works for r — > case because eqs.(j30|) and (|8l| ) do 
not include the mass of the muon explicitly. 

6 Summary 

In this paper, we have calculated the differential cross sections of e + e~ — > t + t~ — > /b/a 
processes, where one of r's decays through LFV processes. Using spin correlations of t + t~, 
we show that the P odd asymmetry of r — > /ry and r — > e7 and P and T asymmetries 
of three body LFV decays of r can be obtained by angular correlations. These P and 
T odd quantities are important to identify a model of new physics responsible for LFV 
processes. 

We have also considered the background suppression of the r — > /ry and r — > e7 search 
by the angular distributions. We see that the analysis of the angular distributions are 
useful for the r + — > fi^y (r~ — > n^l) an d t + — > e^7 (r~ — >■ e^7) searches. 

We can obtain similar information in muon decay experiments if initial muons are 
polarized. Although highly polarized muons are available experimentally, a special setup 
for production and transportation of a muon beam is necessary for actual experiment. 
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The advantage of the r case is that we can extract the information on r spins by looking 
at the decay distribution of the other side of r decay so that we do not need a special 
requirement for experimental setup. 
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A The derivation of the general formulae 

In this section, we derive the eq.(|TJ) from the amplitude in eq.(H). 

By using the completeness relation of the fermion spinors, the amplitude squared is 
deformed to 

\A(p A + m T )^(p B -m T )B\ 2 

2 

E E ^ m (PaAMpa, AihXps, A 2 )u(pb, a 2 )s 

Ai=± A 2 =± 

= E E E E (Mpa,Xi)u(pa,K)A) 

Ai=± A 2 =± A' 1= ± \' 2 =± 

x {u{p A , Ai)tMpb> A 2 Mpb, X' 2 )Yu{p A , A;)) (Bv{Pb, X' 2 )v{Pb, A a )S) , (83) 

where A's are the spin eigenvalues. The spin summation can be performed by using the 
Bouchiat-Michel formulae as follows IP ]: 

\A(f) A + m T )^(i) B -m T )B\ 2 
= a D ~ Tr [(J) a + m T )Y(f> B ~ m T )j v ] a D + 

+ a D - Tr [(j> A + m T )nMB ~ m T )j v ] pf + 
+ pf - Tr ftsfiStA + m T )j^ B - mj-f) oP+ 
+ p D a - Tr [-KfitfA + m T )YlAd>B ~ mj-f] p\ 



where 



a D ~ = ~ {A(t> A + m T )A) , a D+ = i {B(j> B - m T )B) 



(84) 
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Pa - = ~ {A 75 A {fa + m T )A) , p D h + = i {5 75 ^ - m r )i?} , (86) 
where and (s b B ) u are four vectors which satisfy the following equations: 

PA-s a A =p B -s b B = 0, (87) 

-b „a „b raft foo\ 



b A b A — b B b B — 



£(^0$), = + , E(4),(4), = + • (89) 

a=l T 6=1 

The second and third terms in eq.(^4[) vanish because the production parts are antisym- 
metric on fi and v indices while the square of the electromagnetic current from e + e~ 
collision is symmetric on \x and v indices. Explicit calculation gives 

Tr [{fa + mrW^Ufa - m T )Y] = Um T e^ VBp { b B \ + 4im T e^p Ap (s b B ) a , (90) 
Tr [jMfa + m T )Y(fa - m T )r] = Um T e^ PBp (s a A ) a + Aim T e^p Ap {s a A ) rj . (91) 

^ \Ve+lnUe-\ 2 = Tr \t> e +lpi>e-lv\ 



spin 



^Pe+yPe-u + 4p e+u Pe- p - 4^p e + • p e - . (92) 



Using the narrow width approximation, 



1 



q 2 — (m — y)' 



2 



71 



mT 



5(q 2 - m 2 ) , (93) 



the first and last terms in eq. (j84T) give formula (|l|) after the phase space integral. 

B The kinematical functions 

In this section, we list the kinematical functions used in the formulae of branching ratios. 

The functions «i_5, /3i-2, and 71-4 in the r + — > ji + fi + \T and r + — > ji~e + e + decay 
calculations are given as follows. These functions are the same as those used in /i + — > 
e + e + e~ decay fl2"l . x\ and x 2 are given by x\ = 2E\jm T and x 2 = 2E 2 /m T . 

a 1 (x 1 ,x 2 ) = 8(2 — xx — x 2 )(xi + x 2 — 1) , (94) 
a 2 (x 1 ,x 2 ) = 2 {27(1 - xi) + x 2 (l - x 2 )} , (95) 
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f 2x1 - 2x 2 + 1 2x\ - 2x 1 + 1 1 
«.(«,,*) = 8{ \_l + }, (96) 

a 4 (xi, x 2 ) = 32(xi + x 2 - 1) , (97) 

a 5 (xi,x 2 ) = 8(2 - xi - x 2 ) , (98) 
2(x 1 +x 2 )(x? + x2)-6(x 1 +x 2 ) 2 + 12(xi+x 2 )-8 

Pi(xi,x 2 j — (yyj 

2 — xi — x 2 

Q 

p2{Xl,X 2 ) 



(1 - Xi)(l - x 2 )(2 - Xi - x 2 ) 

x { 2(xi + x 2 )(x\ + x\) - A(x 1 + x 2 )(2x 2 + x x x 2 + 2x 2 ) 

+ (19x 2 + 30xix 2 + 19x 2 ) - 12(2xi + 2x 2 -l) > , (100) 



= V(i-^)(i-^)(^ + ^-i)(^-^) (101) 

2 — Xi — x 2 



72(^2) = nJ * 1 ^**-* , , (102) 



/ Xi + x 2 • 


- 1 


y(i-*i)(i 


-£2) 


/ Xi + x 2 ■ 


- 1 


V (i 


-x 2 ) 


/ Xi+ x 2 - 


1 , 



73(^1,^2) = 16 W 7: rrz r (xi + x 2 - l)(x 2 - Xi) , (103) 



7 4 (x 1? x 2 ) = 8y ^ _ x ^^_ x ^ (2-Xi -2:2)^2-^1) • (104) 
The functions Ai_ 7 , -81-4, and Ci_6 in the r + — > yU + e + e~ decay calculation are given 

by 

^ 1 ,. 2 ) = 8(2 -" i - 4 f + 2 ^+ 2 *s, (105) 

1 — Xi 

A 2 (xi, x 2 ) = -8 (xi + x 2 - 1) , (106) 
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A 3 (x 2 ) = -8 (l-x 2 ) , (107) 

m*i) = xi{l : xi) , (los) 



A 5 ( Xl ,x 2 ) = 2 (2 - xi - x 2 ) (xi + x 2 - 1) , (109) 

A 6 (x 2 ) = 2x 2 (l-x 2 ) , (110) 

A 7 (xi,x 2 ) = -8 (4-5xi + x 2 -8x 2 + 4xiX 2 + 4x 2 ,) , (111) 

A 8 (xi,x 2 ) = -4 (l-xi) (xi + 2x 2 -2) , (112) 
#i(xi,x 2 ) 



(1 - xi) (2 - xi - x 2 ) 
x (-6 + 8xi - 3xj + 12x 2 - llxix 2 + 2x 2 x 2 - 8x^ + 4xix 2 + 2xf) , 

(113) 

-(1-m) (2-2xx + x 2 -2x 2 + XiX 2 ) 
B2{Xu X2) ~ 2(2-X!-x 2 ) ' (114) 

2 (l-x 2 ) (2-2x 1 -2x 2 + x 1 x 2 + xl) 

S 3 (xi,x 2 ) = , (115) 

2 - xi - x 2 



5 4 (xi,x 2 ) 



2 — xi — x 2 

x ( - 10 + 16xi - 7x1 + + 22x 2 - 23xix 2 

+5x 2 x 2 - 16x2 + 8x1X2 + 4x 3 2 ) , (116) 

4 (l-xi) (2-4x 1 + x 1 2 -4x 2 + 3xix 2 + 2x 2 2 ) 
S 5 (xi,x 2 ) = 2-X!-x 2 ' (U7) 

(1 - xi) (2-xi -x 2 ) 



c 2fa ,, 2) = ^i±^l^B B EES, (119) 

1 — Xi 
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c 3 ( Xl , X2 ) = - 8 {1 ~ X2) v^-^H 1 -^) EjEtfES (120) 



C 4 (x 1; * 2 ) = {1 ~ X2) ^ + X2 ~^ , (121) 

2 — X\ — X 2 

C 5 (x u x 2 ) = 4 (1 = " 2) ^ = " l} (1 = X2) {Xl +X *- l \ (122) 

2 — xi — x 2 

C e ( Xl , x 2 ) = zWOEg C 1 -^) (*i +**-!) (3-*i-2s 2 ) 

2 — xi — x 2 

C 7 (x 1 , x 2 ) = 8 y/(l - xT) (1 - x 2 ) ( Xl + x 2 - 1) . (124) 
Finally, the functions F, G, and i7 in the r — > /ii/z>y decay calculation are given by 
F = F<°> + rF^ + r 2 F^ , (125) 

G = G (o) +rG W +r i G (2) , (126) 

H = H^+rH^+r 2 H^ , (127) 

where F^ 2 \ G^ 2 \ and #(°M 2 ) are the functions of x(= 2E II /m T ), y(= 2E 1 /m T ), 
d(= 1 + fi^cosz) with ftp = a/1 — Ar/x 2 (r = m^/m 2 ) and z = 7r — 6^ 7 . These functions 
are given by 

(0) -8 (-3 + 2x + 2y) (2x 2 + 2xy + y 2 ) 
F y, (x,y,d) = 

+8x {x 2 (2 + 4y) + y (-3 + y + y 2 ) + x [-3 + y + Ay 2 )} 

-2x 2 y {-Q + y (5 + 2y) + 2x {4 + 3y)} d 

+2x 3 y 2 (2 + y) d 2 , (128) 

(1) 32 (x + y) (-3 + 2x + 2y) 8 {6x 2 + (6 - 5y) y - 2x (4 + y)} 
(X ' y ' rf) = ^ + d 

_ 8x {-A-(-3 + y)y + 3x (l + y)} + Qx 2 y (2 + y) d , (129) 

^ (g>y><Q= -32(-4 + 3x + 3 y) + 48y 

a; d z d 
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G(0) = -8x {Ax 2 + y {-l + 2y)+x (-2 + 6?/)} 

+4x 2 {-2 + 3y + 4y 2 + x (4 + Qy)} - 4x 3 y (2 + y)d, (131) 

JX 32(-l + 2x + 2y) 8x(Qx-y) „ . . , 

GW(x,y,d) = —± - 2 ^ + y d y) -12x 2 (2 + y) , (132) 

— Q6 

G( 2 >(x,</,<i) = -|^, (133) 

^ ^ = -«I> (* + ■>) (-1 + 2, + 2,) 

+4xy {2a; 2 + 2y (1 + y) + x (-1 + 4y)} 

-2x 2 y 2 (-1 + 4x + 2y) d + 2x 3 y 3 d 2 , (134) 

(1) _ 32y (-l + 2s + 2y) _ 8 y (-2 + z + 5y) 

W ' j " xd 2 d 

-4x (3x-2y) y + Qx 2 y 2 d , (135) 

= ^ + ^. (136) 
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